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^invention a pour objet un precede de simulation numerique d'une 
interaction non lineaire entre une source excitatrice et une onde dans un 
resonateur et peut s'appliquer, en particulier, a la synthese numerique, en 
temps reel, d'un phenomene oscillant tel que le son emis par un instrument 
5 de musique fonctionnant en oscillations entretenues, comme un Instrument a 
vent. 

On sait que la simulation numerique d'un phenomene auto-oscillant 
peut s'effectuer par discretisation, dans le domaine temporel, d'equations 
constituant le module du phenomene tel qu'il se traduit physiquement Dans 

10 ce cas, cependant, la synthase sonore, en temps rSel, est difficile a r§aliser 
et c'est pourquoi, depuis quelques ann6es, on a d^veloppe des methodes du 
genre « guide d'ondes numeriques », basees sur un formalisme « traltement 
du signal » de la propagation des ondes dans les deux directions du corps de 
{'instrument Habituellement. de telles methodes utillsent la decomposition de 

IS la vibration au sein du resonateur sous forme d'une somme des ondes aller^ 
retour se propageant. respectivement, de la source ^ rextr§mit6 dui 
resonateur et de cette extremite a la source. Une telle methode ne permet, 
pas d'exprimer simplement le couplage non lineaire qui existe entre la source*, 
excitatrice et le resonateur de Tinstrument et limite la parametrisationv 

20 physique des algorithmes de synthase. De plus, une telle m§thode limite le 
choix de la geometrie du corps du resonateur a des portions cylindriques de 
tubes et rend difficile la prise en compte de la conicite du resonateur, par 
exemple pour des instruments du type saxophone, 

Uinvention a pour objet de resoudre ces inconvenients et de lever ces 

25 limitations grace a une nouvelte methode de simulation et de synthese en 
temps reel d'un phenomene oscillant, applicable specialement mals non 
limitatlvement aux instruments a vent. En particulier, Tinvention a pour objet 
un proc§de de simulation permettant de prendre en compte les processus 
physiques r6gissant le fonctionnement d'un instrument reel et dont la mise en 

30 oeuvre numerique peut etre particulierement simple. 

En outre» a partir d'une methode de base applicable a des Instruments 
du type clarinette, a resonateur cylindrique, Tinventlon peut §tre adapt§e ^ la 
simulation d'autres types d'instruments. 
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Par allleurs. I'invention ne se limite pas d ia simulation d'instruments 
de musique mais peut s'appliquer, d'une fagon g6nerale d la synthese 
num^rique, en temps reel de phenomenes oscillant de toutes sortes. 

D'une fagon generate, I'invention a done pour objet la simulation d'une 
5 Interaction non lineaire entre une source excitatrice et une onde dans un 
resonateur, au moyen d'outils de calcul de signaux numeriques, soit ^ partir 
d'equations dont la solution correspond a la manifestation physique d'un 
phenomene a simuler,, soit a partir d'une mesure de cette manifestation 
physique. 

10 Conform6ment a I'invention, le phenomene a simuler se traduisant, a 

chaque instant et en un point donn§ du resonateur. par une relation 
d'lmp6dance entre deux variables representatives de I'effet et de la cause 
dudit phenomene, on transcrit directement I'^quation de I'lmpedance sous 
fomie d'un module faisant appel a des guides d'onde pour r6allser une 

15 interaction non lineaire entre les deux variables de la relation d'impedance. 

A cet effet, le modele comporte, d'une part, au moins une partle 
lin§aire representant directement rimp6dance dite d'entr6e du resonateur, 
c'est-^-dire rimp§dance au point oCi se manifeste I'interaction non lineaire et, 
d'autre part, une partie non lineaire modeilsant le r6le de la source excitatrice 

20 du phenomene a simuler. 

En partlculier, pour la synthase num§rique, en temps reel, d'un 
phenomene oscillant, I'invention permet, d partir d'un systeme d'equations 
entre au moins deux variables representatives du comportement du 
resonateur, d'etablir une expression de I'impedance d'entree du resonateur 

25 en terme de guide d'ondes lineaires, sans decomposition en ondes aller- 
retour, de fagon a realiser au moins une partie lineaire du modeie qui peut 
etre coupiee S une boucle non lineaire faisant intervenir revolution de la non- 
linearite telle qu'elle s'exprime entre les deux variables de la relation 
d'impedance d'entree du resonateur. 

30 De fagon particulierement avantageuse. cette partie lineaire du 

modeie constitue ia somme de deux guides d'ondes eiementaires realisant 
une fonction de transfert entre les deux variables de la relation d'impedance. 

Selon une autre caracteristlque particulierement avantageuse, le 
modeie est pilote par au moins deux parametres representatlfs de 

35 I'interaction physique non lineaire entre la source et le resonateur, au moyen 



d'une boucle reliant la sortie a Fentree de la partie lineaire et comprenant une 

fonction non lineaire jouant le role de source excitatrice pour le resonateur. 

Comme on Fa indique, Tinvention s'applique, en particulier, a la 

synthese en temps reel des sons produits par un instrument a vent. Dans ce 
5 cas, les deux variables de la relation dimpedance sont la pression et le debit 

acoustiques d Tentree du resonateur. 

Dans le cas d'un resonateur cylindrlque ayant une extremite ouverte, il 

est particulierement avantageux de r§aliser la partie lineaire du module de 

transcription nunnerique de r§quation de rinnp6dance sous forme d'une 
10 somme de deux guides d'ondes elementaires ayant pour source excitatrice le 

debit a Tentr^e du resonateur et r^alisant la fonction de transfert : 

Z { ) ^ ^^^^^ = ^ exp(-2ik((D)L) 

^^^^ Ue(a>) ~ 1 + exp(-2ik(©)L) 1 + exp(-2ik(a>)L) 

15 dans laquelle : ^ 

' ' «^ 

- 00 est la pulsation de Tonde 

- Ze(a)) est Timpedance d'entree du resonateur, ^ 

- Pe(co) et Ue(o)) sont les transformees de Fourier des valeurs;^, 
adimensionnees de la pression et du debit a I'entree du^ 

20 resonateur, 

- k((D) est une fonction de la pulsation de Tonde qui d6pend du 
phenomene S slmuler, 

- L est la longueur du resonateur. 

Selon une autre caracteristique, chacun des deux guides d'ondes fait 
25 intervenir un flltre ayant pour fonction de transfert : 

-F(a))^ =-exp(-2ik((]o)L) 

et representant un trajet aller-retour d'une onde avec changement de 
30 signe ^ Textremite ouverte du resonateur, chaque guide d'onde 
correspondant a un terme de r§quatlon de Timpedance. 

Un tel module peut avantageusement §tre pilote par la longueur du 
resonateur et au moins deux parametres representatifs de Tinteraction 
physique non lineaire entre la pression et le debit a Tentree du resonateur, au 



moyen d'une boucle reliant la sortie d I'entr^e de la partle llneaire et 
comprenant une fonction non Iln6alre jouant le rdle de source excitatrice pour 
le resonateur. 

En particulier, pour la synthase en temps r§el des sons S simuler, on 
5 realise une formulation, dans le domaine temporel, de la reponse 
impulsionnelle du resonateur, en approximant les pertes representees par le 
filtre au moyen d'un filtre num^rique approche. 

L'inventlon couvre d'autres caracteristiques essentlelles mentionnees 
dans les revendlcations et portant, en particulier, sur les equations utilisees 
10 par I'outll de calcul des signaux numeriques et qui conduisent a des modeles 
^ guides d'ondes dependant du phenomene d simuler. 

En effet, selon une caracteristique essentielle de I'invention, la 
methode proposee pour la simulation d'un phenomene simple comme la 
propagation d'une onde dans un resonateur cylindrique, peut etre adaptee de 
15 multiples fagons pour la simulation de phenomdnes plus compiexes et, en 
particulier, de divers types d'instruments. 

Dans la description qui va suivre, on va done exposer en detail la 
methode de simulation, les equations utilisees et le modele d mettre en 
oeuvre pour la synthese du son d'un instalment d resonateur acoustique 
20 cylindrique a anche, du type clarinette. ainsi que certalnes adaptations pour 
la simulation d'autres types d 'instruments. 

La figure 1 donne deux diagrammes representant respectlvement, ^ 
gauche la fonction de transfert, en Hertz, d'un modele d'anche ^ un mode et 
a droite, la reponse impulsionnelle en fonction des echantillons, avec une 
25 frequence d'echantillonnage de 44 100 Hertz. 

La figure 2 est un schema de calcul par combinaison de guides 
d'ondes, representant I'impedance d'entree d'un resonateur cylindrique. 

La figure 3 donne deux diagrammes representant respectlvement, 
pour un resonateur cylindrique, en haut I'impedance d'entree en fonction de 
30 la frequence indiquee en Hertz et. en bas, la reponse impulsionnelle en 
fonction du temps, en secondes. 

La figure 4 est un schema de calcul d'un modeie de simulation d'un 
instrument a anche ^ resonateur cylindrique. 

La figure 5 donne deux diagrammes analogues i la figure 3, 
35 representant respectlvement, pour un modeie de resonateur calcuie selon 



invention, en haut rimpedance d'entree approchee et en bas la reponse 
impulslonnelle approchee. 

La figure 6 donne deux diagrammes analogues a la figure 1 • 
repr§sentant respectivement. pour un modele d'anche calcule selon 
rinvention. a gauche la fonctlon de transfert et a drolte la r6ponse 
impulslonnelle. 

La figure 7a montre les variations, en fonctlon du temps indique en 
secondes. de la pression acoustique Interne au niveau de I'embouchure d'un 
resonateur cylindrique. Les figures 7b et 7c sont des agrandlssements des 
transitoires d'attaque et d'extinction. 

La figure 8 donne deux diagrammes representant respectivement d 
gauche la fonction de transfert et ^ drolte la reponse Impulslonnelle, pour un 
modele d'anche multimode calcule selon I'Invention. 

La figure 9 donne deux diagrammes representant le spectre de la 
pression acoustique externe, respectivement en haut pour.une anche S upt 
seul mode et en bas pour une anche ^ modes multiples. t 

La figure 10 est un schema de calcul representant I'impedance d'unt 
r§sonateur cylindrique avec impedance terminate. f 

La figure 1 1 est un schema de calcul representant rimp§dance d'un? 
r§sonateur conique. 

La figure 12 est un schema de calcul representant I'impedance d'un"" 
resonateur simplifie pour cuivres. 

La figure 13 est un schema general de calcul representant 
I'impedance d'une combinaison en paralieie de r^sonateurs cylindriques. 

Comme on le salt, les phenomenes physiques mis en jeu iors de la 
production du son de la clarinette. s'expriment d'une part, par une equation 
de propagation lineaire des ondes dans le tuyau avec perte. et d'autre part 
par une equation non lineaire reliant le debit avec la pression et le 
deplacement de I'anche au niveau de I'embouchure de Tinstrument 

Un modele de simulation du son comprend done une partie Ilneaire du 
modele correspondant au resonateur de I'instrument qui. dans le cas de la 
clarinette est constitue d'un tube cylindrique. 

Pour cette geometrie. en supposant que le rayon du tube est grand 
devant les epaisseurs des couches limites. la pression acoustique § 
rinterieur du tube est regie par une equation de la forme : 



a'p(x,t) 1 a^p(x,t) ^a|p(x^_Q 

oCi « = ]^^(^V^+[~~^]"N^]-' ^ '® "^^y®" c'est-a-dire 
7mm dans le cas de la clarinette. Les valeurs des constantes physiques, en 
unites mKs. sont : c=340. lv=4.10-^ lt=5.6.10-«, •^ = 1-4. 

En recherchant les solutions du type exp (i{a)t-k(Q))x)). o etant la 
pulsation de I'onde, on peut 6crire : 

En rempla9antVl+x par la valeur approchee 1 + |- lorsque x est 

petit, I'expression approchee classlque de k(<D), que nous utillserons par la 
suite, devlent : 

k(co) = ^-i^ac«>^ (2) 
c ^ 

On sait que, si Ton consldere un tuyau de longueur infinie que I'on suppose 
excite en x=0 at t=0 par une impulsion Dirac 5(x)5(t), en tout point x>0, la 
pression acoustique propag6e a partir de cette source s'6crit sous la forme 
d'une somme continue de toutes les ondes susceptibles de se propager dans 
ie tuyau : 

p(x,t) = /exp(-ik(a))x)exp(i(Dt)do> 

qui apparatt comme la Transform6e de Fourier Inverse de la valeur 
exp(-ik(<D)). 

La fonction de transfert d'un tuyau de longueur L, qui constitue ie filtre 
du module guide d'onde representant la propagation, la dissipation et la 
dispersion est done : 
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F(<o) = exp(-ik(©)L) = exp|^- ^ ac^ x exp - '^f «- + 1 J (3) 

La dissipation representee par le module de F(o)) et la dispersion 
representee par la phase de F(g)) sont done proporlionnelles § -Ja , tandis 

que le retard de propagation est ddnne par -. La longueur du tuyau sera 

done le parametre de controle de la hauteur et son rayon contrSIe les pertes. 

On sait, d'autre part, que les transfomn6es de Fourier des pressions et 
debits adimensionnes ^ I'entree (Pe(ci)), UeC©)) et en sortie (Ps((o), Us(co)). du 
resonateur, sont lies par le systeme d'equations : 

Pe(a>) = cos(k(co)L)Ps(o))+ / sln(k(<D)L)Us(c(>) 

Ue(co) = /sin(k(©)L)Ps(oi>)+ cos(k(<D)L)U8(co) 



15 - De fagon classique. afin de modeiiser ia pression acoustique interne;^ 
le rayonnement peut etre neglige. L'extremite ouverte de I'instrument estj 
done parfaitement refiechissante, ce qui entraTne que Ps((o)=0. Ceci permet: 
d'exprimer la relation entre pression et debit a I'entree du resonateur : 



Pe(<D) = /■ tan(k((o)L)Ue(co) = Ze (a))Ue(<x)) (4) 



Gill Ze(o))= /■ tan(k(<D)L) est I'impedance d'entree nomnalisee du 
resonateur. 

Dans le cas d'un modele classique d'anche ou de levres a un mode, le 
25 deplacement adimensionne x(t) de I'anehe par rapport d son point d'equilibre, 
et la pression acoustique Pe(t) qui le produit sont lies par I'equation : 

-l^ + i^^x(t) = Mt) (5) 
00? dt=^. dt 



30 



dans laquelle a)r=27c/r correspond a la frequence de resonance /r, par 
exemple 2500 Hz et qr est le facteur de qualite de Tanche. par e>(emple 0,2. 



En ecrivant. I'equation (5) dans le domaine de Fourier, on obtient la 
fonction de transfer! de I'anche : 

X((o) _ of 
Pe(co) cof - o>2 + icoq,(D^ 

dont la r§ponse impulsionnelle est donnee par : 
^<*> = ^/^qr^ ®''P(- I «*r^rt] sin[^l ^^c.t j (7) 

Comme on le verra plus loin, il est important de noter que x(0)=0. 

Par ailleurs, dans le cas d'un instrument a anche de type clarinette ou 
S embouchure de type trompette. la pression acoustique pe(t) et le d^blt 
acoustique Ue(t) (adimensionnds) S {'entree du resonateur sont relics de 
fagon non-lineaire par I'equation : 

Ue(t) = (1 - sign(y - x(t) - 1))sign{Y - Pe(t)K(1 - y + x(t)y|y - p,{t)| (8) 



Dans le cas d'un Instrument ^ anche. le paramdtre ^ est 
caract^ristique de I'embouchure et prend en compte la position des levres et 
le rapport de section entre le bee et le resonateur. Des valeurs classlques 
peuvent etre de 0,2 a 0.6. Le param^tre y est !e rapport entre la pression ^ 
rinterieur de la bouche d'un instrumentiste et la pression de plaquage 
statique de I'anche. Pour un tuyau sans perte, ii va de i pour la mise en 

vibration a 1 pour la position d'anche battante. 

Les paramdtres ^ et y sont done deux parametres de jeu Importants 
dans la mesure oO ils repr6sentent, respeetivement, la maniere dont 
rinstrumentlste pince I'anche et la pression du souffle dans I'instrument. 

En comblnant {'equation de deplacement de I'anche, la relation 
d'Impedance et la caracteristique non lineaire. 11 apparaTt que la pression 
acoustique au niveau de I'embouchure est contrdlee par le systeme 
d'equations suivant : 
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1 d^x(t) ^ q, dx(t) 
cof dt' 00^ dt 



+ x(t) = p,(t) 



(9) 



Pe(a>) = i tan 



— -^^ac(D''Ou,(a>) 
c 2 ^ 



(10) 



5 Ue(t) = -(1 - sign(y - x(t) - 1))slgn(y - Pe(t)K(1 - Y + x(t))^|y - p, (t)| (11) 
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.15 



20 



La figure 1 donne deux diagrammes indiquant respectivement, d 
gauche la fonction de transfer! et a droite, ia r§ponse impulsionnelle du 
modele d'anche, pour une frequence de resonance /p= 2 500Hz et un facteur 
de quality qr=0,2. 

L'id§e de I'invention est done de trouver une formulation dans le. 
domaine temporel de la relation d'imp^dance permettant de resoudre ce ■ 
systeme de trois Equations, en mod^lisant la relation d'imp§dance en termes ' 
de modeles de guides d'onde elementaires. 

Pour modeliser Timpedance d'entree du resonateur en terme devi 
modeles de guides d'onde elementaires, on ecrit la transformee de Fourier ' 
de I'impedance Ze(co) sous la forme : 

ZJ<o) = itan(k(a))L) = i S'"(k(fi>)'-) ^ exp(il<(a))L) - exp(-ik(o))L) 

cos(k(ffl)L) exp(ik((D)L) + exp(-ik(o))L) 

Cette expression peut s'6crire sous ia forme : 



2,(0)) = 



Pe((D) ^ 1 exp(-2ik((a)L) 

Ue(Q)) 1 + exp(-2ik((D)L) 1 + exp(-2ik(a))L) 



(12) 



25 La figure 2 represente un modele de calcul par combinaison de guides 

d'ondes directement issu de cette derniere Equation et dont la fonction de 
transfert est I'impedance d'entree du resonateur. II est constitu6 d'une 
somme de deux guides d'ondes elementaires. L'el^ment superieur 
correspond au premier terme de I'^quation (12) tandis que I'element inferieur 

30 correspond au deuxieme. Le filtre dont la fonction de transfert est 



- F((a)^ = - exp(-2ik(cD)L) represente un trajet aller-retour. avec changement 
de signe de la pression acoustique a rextremite ouverte. 

A titre d'exemple, la figure 3 donne deux diagrammes representant 
respectivement, pour un tuyau de longueur L = 0.5m et de rayon R=7mm, en 
5 haut la variation de ('impedance d'entree du r^sonateur en fonction de la 
frequence et, en bas, la r6ponse impulsionnelle du modele a guides d'onde 
correspondant, calculee par transformee de Fourier inverse de rimp^dance. 

Le syst^me des trois Equations physiques couplees (9), (10), (11) 
permet d'introduire la non-linearlte sous fonne d'une boucle reliant la sortie 
10 pe du rSsonateur d I'entrde Ue. 

La figure 4 donne un schema de calcul Equivalent permettant, pour la 
simulation d'un instrument S anche ou embouchure, de coupler non 
lin^airement le d§placement de I'anche et la pression acoustique avec le 
debit acoustique a I'entree du r6sonateur, en calculant, d chaque instant 
15 echantillonn§, la pression acoustique interne au niveau de I'embouchure. 

Le moddle est entierement pilots par la longueur L du rEsonateur et 
au moins deux param6tres ^ et y repr6sentatifs de interaction physique non 
lineaire entre la. source et le resonateur, au moyen d'une boucle reliant 
I'entree a la sortie de la partie lineaire et comprenant une fonction non- 
20 lineaire jouant le role de source excitatrice pour le r§sonateur. Comme le 
montre la Figure 4, la partie lineaire reprend le schema de la Figure 2 et la 
fonction non-lineaire f est controlee par les deux param^tres ^ et 7 
pemnettant de simuler le jeu d'un instrumentiste, et a comme parametres 
d'entree, dans le cas d'une clarinette, la pression a I'embouchure et le 
25 deplacement x(t) de I'anche par rapport ^ son point d'equilibre, calcule en 
fonction de la pression S I'embouchure. par un modele d'anche (m). 

Aucun signal d'entree n'est n^cessaire, dans la mesure oCj y est 
difBctement proportionnel d la pression dans la bouche de I'instrumentiste. 
C'est done la non-lin§arit§ elle-mSme, et son Evolution Imposee par 
30 I'instrumentiste, qui joue le r6le de source excitatrice, conformEment au 
module physique. 

Pour la synthese. en temps reel, des sons d simuler, le module 
necessite un echantillonnage numerlque et, pour cela, on realise une 
formulation, dans le domaine temporel, de la rEponse impulsionnelle du 
35 rEsonateur. correspondant a la transformEe de Fourier inverse de 



I'impedance. Cette formulation dans le domaine tenriporel permet de calculer 
la pression pe(t) a I'embouchure en fonction du d^bit Ue(t) mais il est 
necessaire, pour cela, d'approximer les pertes representees par le filtre F((o) 
au moyen d'un filtre numerique approche. On va done utiliser une 

5 approximation de la fonction de transfert du filtre, - F((a)2 = - exp(-2ik(o)L) , 
dont les coefficients sont determines d partir de variables physiques telles 
que la longueur du r§sonateur et son rayon, de fagon a pouvoir apporter les 
modifications n§cessaires en fonction de la geom§trie du resonateur. A cet 
effet, on exprime analytiquement les coefficients du filtre num6rlque comme 

10 des fonctlons de param§tres physiques. 

En pratique, il est partlcull^rement avantageux d'utillser un filtre § un 
p6te en exprimant la fonction de transfert approch6e sous la forme : 

a, . b- exp(-2itDD) 



15 



1 - a^ exp(-itiy) 



dans laquelle m = j-, /e etant la frequence d'6chantillonnage, eife 

L ^ 
D = fe - est le retard pur correspondant i un trajet alter ou retour des ondes^ 

dans le resonateur. '"^ 
Les parametres bo et ai sont exprim§s en fonction des parametres 
physiques de telte sorte que |F(<»)2|^ = |F(tij)|^ pour deux valeurs donn6es de 

20 CO, 

La premiere valeur ©i retenue est celle de la frequence fondamentale 
de jeu. Ceci perniet d'assurer un temps de decroissance de la frequence 
fondamentale de la reponse impulsionnelle du modele a guides d'ondes 
utilisant le filtre approch§. Identlque ^ celul du module k guides d'orides 
25 utilisant le filtre exact. 

La deuxieme yaleur ©2 retenue est celle d'un harmonique choisi de 
fagon d obtenir une decroissance globale identique des r§ponses 
impulslonnelles des guides d'onde, respectivement, exact et approche. 

Le choix de cette deuxieme valeur (02 est done plus libre. Elle 
30 correspond, par exemple, au second pic de resonance dans le cas de la 
clarinette mais, dans certains cas, comme on le ven-a plus loin dans le cas 
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de la trompette, il peut Stre preferable de choisir un harmonique de rang 
superieur. 

Le systeme d'equations a resoudre est donne par : 
|F(a),)f (1 + a? - 2a, cos{w,)) = b§ 



|F(a>2)f (1 + a,2 - 2a, 003(^2)) = b§ 
dans lequel |F(co)^|^ = exp - Zac^j^L^. 

En negligeant la dispersion introduite par la partie non-liniaire de la 

4-1) 

10 phase de F((d), les frequences des harmoniques sont: tn,- = — -oil 

Lfe 

2k-1 est le rang de I'harmonique. 
En notant : 

c, = cos(t0,),C2 = cos(to2).Fi = |F(o),f |^F2 = |F((D,)''f . A, = FjC^Ag = F^c^, 
les coefficients a^ et bo sont donnes par : 
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a, = A1 - A2 - V(A, - A,)^ - (F, - F,)^ ^^^^ 

F, - F2 



. V2F,F2(c, - C3)(A, -A,- .liKT- A^f - (F, - F^)') 

Do = (15) 



20 A partir de I'expression de I'impedance d 'entree du resonateur 



Ze(o)) = exp(-2ik(co)L) 

1 + exp(-2ik(o))L) 1 + exp(-2ik((D)L) ' 



et en notant z=exp(im), on tire directement 



Pe(z) 
Ue(z) 



( 



,-2D 



-1 



1 + 



-2D 



-1 



1 + 



1 - a,z 



-2D 



1 - a^z 



-1 



a^z 



-1 



+ boZ 



-2D 



d'oO Ton tire l'6quation aux differences : 

5 

Pe(n) = U3(n) - a^UeCn 1) - boU^Cn - 2D) + a,Pe(n ~ 1) - boPe{n - 2D) (16) 
La figure 5 donne deux diagrammes montrant respectivement, en haut 
la variation de Tinnpedance d'entree du resonateur approche en fonction de la 
frequence et, en bas, la reponse impulsionnelle du modele a guides d'onde 
10 correspondant, calculee a partir de Tequation aux differences pour un tuyau 
cylindrique de longueur L=0.5m et de rayon R=7mm. 

On constate que le nnodeie ainsi etabli donne une reponse 
impulsionnelle tres proche de celle du resonateur, representee sur la figure 
1. -n- 
15 De meme que pour le filtre reprSsentant les partes, la relation entre ja 

pression acoustlque et le dSplacement de Tanche doit §tre discretisee dans 
le domaine temporel. Or, la reponse impulsionnelle de I'anche est une 
sInusoTde amortie de fagon exponentielle qui satisfait a la condition x(0)=0, 
comme on Ta indique plus haut. II est done possible de construire un filtre 

20 numerique pour lequel le deplacement de Tanche a Tinstant t^ = — est 

fe 

n — 1 

fonction de la pression acoustique d I'instant t^^^ = et non pas tn. Ceci 

fe 

permet de respecter la propriete x(0)=0 du systeme continu lorsque Tanche 

est soumise a une excitation Dirac. Pour satisfaire cette condition, au lieu 

d'utiliser la transformation bilineaire pour approximer les termes ico et -co^, on 

fe 

25 utilise, selon Tinventlon, les expressions io^— (z~z~^) et 

-co^«fe^(z - 2 -h z'''), qui correspondent a une differentiation numerique 

centree exacte au second ordre. 

Avec ces approximations, la fonction numerique de transfert de 
Tanche est donnee par : 
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X(Z) _ fflf 

P^(^> CD? + U-iz - 2 + z-^) + I (z - z-^)q,a>. 















2®r 




2 


1 2®, Cof; 
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15 



20 



25 



d'oa I'on tire I'^quation aux differences : 

x(n) =0xpe(n)+biaPe(n-1)+aiaX(n-1)+a2aX(n-2) 

dans laquelle les coefficients bia, a^a et aza sont d6finis par 



G>r 2(0, 



1 



2f 2 

7 

(Or 



(17) 



*0a 



a 



a2a= 



Oa 



a„ 



'Oa 



La figure 6 donne, pour un tel modele d'anclie approcii6, deux 
diagrammes representant, ^ gauche la fonction de transfert et a droite la 
reponse impuisionnelle, avec une frequence d'echantillonnage fe=44 100 
Hertz, les valeurs des parametres etant fr=2500Hz et qr=0,2. 

On voit que les diagrammes obtenus sont tres proclies de ceux de la 
figure 1. 

Ceci etant etabii, on va maintenant exposer une methode de 
resolution explicite, selon i'invention, pour !e couplage des Equations aux 
differences avec ia caract§rlstique non lineaire. 

Les formulations 6chantiIlonnees des reponses impulsionnelles du 
deplacement de I'anche et de rimp6dance permettent, en effet. d'ecrire 
r^quivalent §chantillonne du systeme d'equation (9. 10, 11) indique plus 
haut, sous la forme : 



x(n)=biape(n-1)+aiaX(n-1)+a2aX(n-2) 



(18) 



pe(n)=Ue(n)-aiUe(n-1)-boUe(n-2D)+aiPe(n-1)-boPe(n-2D) (19) 

"e(n) = ^ (1 - sign(y - x(n) - 1))sign(y - p,(n)K(1 - y + x(n))^|y - p,(n)| (20). 

Ce syst6me d'equations est implicite, du fait que (e calcul de pe(n) par 
I'equation d'imp6dance necessite la connaissance de Ue(n) et que cette 
valeur est elle-meme obtenue S partir de I'equation non lineaire et n§cesslte 
la connaissance de pe(n). Cependant. comme on I'a lndlqu§ plus haut. le 
calcul de x(n) ne necessite pas la connaissance de pe(n) mais celle 'de 
Pe(n-1) qui est connu a I'instant n. 

Ceci permet. selon I'inventlon, de resoudre simplement et exactement 
le systeme couple. Pour cela, les tennes des Equations 19 et 20 ci-dessus, 
qui ne dependent pas de I'echantillon de temps n, peuvent, en effet, §tre 
regroupes dans les expressions : 

15 V=-aiUe(n-1)-boUe(n-2D)+aiPe(n-1)-boPe(n-2D) 

W = 1 (1 - sign(y - x(n) - 1)) x ^(1 - y + x(n)) " 1 
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ce qui permet d'ecrire : 

Pe (n) =Ue(n)+V 



Pour generaliser la methode, 11 est int^ressant d'associer Ue (n) a un 
coefficient bco = 1 dans le cas d'un tuyau cylindrique. 

25 Les deux equations 1 9 et 20 ci-dessus peuvent alors §tre ecrltes sous 

la forme : 

Pe(n)=boCoUe(n)+V 
30 u^(n) = Wsign(Y - Pe(n))V|Y - Pe(n)| 

Etant donne que le ternie 1(1 - sign(y -:- x(n) - 1)). annule W lorsque 
(1-)H-x(n)) est n6gatif. W reste toujoure posltlf. Si J'on considere 



successivement les deux cas : y-pe(n)>0 et 7-Pe(n)<0 correspondant 
respectivement aux cas Ue(n)>0 et Ue(n)<0, Ue(n) peut s'exprimer exactement 
et sans faire intervenir Tinconnue Pe(n), sous la forme : 

5 Ue(n) = ^ slgn(Y - V)(-bCoW2 + \N^{bc^\N)^ + 4\y - V| 

De ia sorte, ie calcul de la pression acoustique et du debit au niveau 
de Tembouchure a un instant echantillonnS n, peut s'effectuer en utilisant de 
fagon sSquentielle les equations suivantes : 

. 10 x(n)=biape(n-1)+aiaX(n-1)+a2aX(n-2) (21) 

V=-aiUe(n-1)-boUe(n-2D)+aiPe(n-1)-boPe(n-2D) (22) 

W =: 1 (1 - sign(Y - x(n) ^ - Y + x(n)) (23) 

15 

Ue(n) = I sign(Y - V)(-bCoW2 + W^(bCoW)^ + 4|y - V| (24) 

Pe(n)=bcoUe(n)+V (25) 

20 Ainsi, {'invention permet de resoudre dans Ie domaine tempore! Ie 

systeme d'equations regissant la modelisation physique de Tinstrument, a 
partir d'une formulation echantillonnee equivalente de la reponse 
impulsionneiie du deplacement de ranche, de la relation dimpedance et de 
la caracteristique non lineaire, qui se traduit par Ie systeme d'equations 18, 

25 19. 20, dans lequel : 

- Tequation (18) est une transcription numerique du modele (m) de la 
figure 4, 

- Tequation (19) est une transcription num6rique du modele 
d'impedance de la figure 2, 

30 - Tequatlon (20) est une transcription num6rique de la 

caracteristique non lineaire reliant Ie d§placement de Tanche et la 
pression acoustique avec Ie debit acoustique. 



En regroupant les termes qui ne dependent pas de I'^chantillon 
temporel n, le precede selon rinventlon permet. en effet. de determiner le 
debit et la pression a I'entree du resonateur par un calcul sequentiel des 
equations 21 a 25, et de resoudre, dans le domaine temporel. le systdme 
d'equations 9. 10. 11 regissant la modelisation physique d'un Instrument a 
anche du type clarlnette. afin de synthetiser les sons produits par un tel 
instrument. 

A titre. d'exempie, la figure 7 montre la variation de la pression 
acoustique Interne au niveau de I'embouchure. calculee par un tel module 
non Im^aire falsant appel S des guides d'ondes. pour un tuyau de longueur 
L-0.5m et de rayon R=7mm. les valeurs des parametres etant y=0 4 C=0 4 
fr=2205 Hz. qr=0.3. . 

On observe trois phases : le. transltoire d'attaque correspondant a une 
augmentation brutale de y et ^. le regime permanent pendant lequei y et r 
diminuent progressivement. de fagon linealre, jusqu'au seull d'oscillation et 
le transltoire d'extinction. * : 

En pratique, la mise en oeuvre num6rique d'un tel module a guides 
d'ondes non lineaire. peut se faire en utilisant des elements setrouvant de^ 
le commerce. Par exemple, on peut effectuer une mise en oeuvre num^riqUe 
en langage C sous forme d'un objet externe « clarlnette » pour 
lenvironnement connu sous la marque Max-MSP, pilote 3 partir de 
commandes MIDI foumies par un controleur Yamaha WX5®. Ce contr6leur 
mesure la pression des l§vres sur I'anche. qui controle le parametre ^ et la 
press,on du souffle, qui controle le parametre y. Ces informations revues au 
format MIDI (done entre 0 et 127) sont renormalis§es pour correspondre ^ 
i'^chelle des parametres physiques. L'accord du guide d'onde est realise a 
partir de I'informatlon MID! pitch contr6lee partir du doigte qui determine la 
longueur L du tuyau. 

Cependant, comme dans un Instrument reel, la hauteur change en 
fonchon des parametres physiques tels que y. C co, et q,. Or. en realite. un 
instrument de musique n'est pas parfaitement accorde pour tous les doigt6s 
L'utilisation d'un filtre passe-tout pour la mise en oeuvre de la partie 
fractionnaire du retard D n'est done pas indispensable. 

En pratique, il apparalt que les sensations de jeu d'un tel instrument 
virtuel sont tout e fait comparables d celles d'un Instrument reel. 



Cependant, le procede qui vient d'etre decrit pour la simulation du son 
d'un instrument a anclie du type clarinette, peut encore etre perfectionne. 

On salt, en effet, que la pression acoustique au niveau de 
I'embouchure n'est pas la variable representative du son pergu. il est done 
5 interessant de calculer la pression externa qui, pour un tuyau cylindrique, 
peut s'exprimer comme la d§rivee en fonction du temps, du debit sorlant : 

Pexi(t)= ^^^^^^ > En negligeant encore le rayonnement, ce qui entraTne Ps(t)=0, 
dt 

il vient : 

10 Pe(<j£)) = / sin(k(a))L))Us(co) 

Ue(OE)) = COS(k(0))L)Us((D) 

d'ou Ton tire : 

Us(<d) =exp(Mk(co)L)(Pe(a))+Ue(a))) 

15 

Du point de vue perception, le terme exp(-ik(co)L) est negligeable, 
L'expression ci-dessus peut done etre simplifiee et devient : 

Pext(t)=;^(Pe(t) + Ue(t)) (26) 

at 

20 

Ainsi, d*un point de vue numerique, le calcul, a chaque instant 
echantillonne (n), de la pression externe Pext(n), se reduit a une simple 
difference entre la somme de la pression interne et du debit, entre i'instant 
(n) et rinstant (n-1). 

25 ' D'autre part, des mesures physiques ont montre que les vibrations 
d'une anche etaient plus complexes qu'une simple sinusoTde amortie. II est 
done apparu que le modele simple decrit plus haut pouvait etre perfectionne 
de fagon a am6liorer la qualite du son produit, tel quil est pergu. Pour cela, 
on considere un modele tres simplifie d'anche sous forme d'une corde 

30 encastree libre. Ainsi, le module d'impedance d'un tuyau cylindrique decrit 
plus haut, qui g§nere principalement des harmoniques impaires, sert de base 
pour la realisation d'un module d'anche a mode multiple satisfaisant la meme 



condition x(0)=0, ce qui permet de conserver le. scliema de resolution 
nurti^rique des Equations 21 a 25 indiquees plus liaut. 

Etant donne que la valeur de I'impedance est rieelle pour tous les pics 
de resonance, sa reponse impulsionnelle est une somme de fonctions 
cbsinus. D'autre part, on a vu que la reponse impulsionnelle d'un modele 
d'anclie a modef unique est une fonctlon sinusoYdale. On va done utillser, 
dans la definition du module, la fonction de transfert entre une sInusoTde 
amortie et une cosinusoTde amortie. En conservant les m§rries notations, la 
transformee de Fourier de x(t)=--^=L= exp^- ^ co^q^tj sin^-^ ^4 - qr^co^tj 



10 est donnee par : X(co) = 



De 

y(t)= exp^- 1 co,q,tj cos|^^ ^4 - q>,tj est donn6e par : 



mSme, 



la 



ioq^oor 
transformee 



de 



Fourier 



de 



IS 



Y(a>) = 



<D,(<o,qr + 2\(o) 



74 - qf (oof - + mqrO>r) 



La fonction de transfert entre y(t) et x(t) est alors 



X(co) ^ V4-qf(a, 
Y(ctt) (s>,qf + 2\(o 

On peut ainsi 6crire le module de transmittance sous la forme : 
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X(a)) 
Pe<o)) 



1 - B exp - 1 — 



( 



-a),q, +1(0 



1 + B exp 



(_, 



Pour determiner les trois inconnues ©a, C, B, on s'Impose trois 
conditions. La premiere est de conserver la frequence du premier pic, qui est 
le pic maximal. On choisit pour cela, a)a= cor en supposant que le d^calage de 
frequence resultant du facteur de quality qr est n^gligeable. 



La seconde condition, satisfalte par le module d'anche S mode unique 
consiste a imposer une valeur unitaire du module de transmrttance pour a>-0 
de fagon a maintenir X(a)) « Pe(a)) aux basses frequences. 

La troisidme condition est une valeur impos6e de — pour le module 

de la transmittance a la frequence cor. afin de conserver la hauteur du pic du 
modele d'anche a mode unique. 

Grace a ces trois conditions, la fonction de transfert ainsi realis§e 
reproduit les caracteristiques principales du modele d'anche ^ mode unique 

Les deux premieres conditions conduisent au systeme d'equations : 

2C(1 + B) 1 
«>rVqr +4(1-B) ~ qT 

2C(1-B) _^ 
Q>rqr(1 + B) 

en posant : A, = . J A2=^<o,q, et A3=A2Aiqr, les coefficients 
B et C r^solvant le systeme sont donn§s par : 
1 + A3 -27a7 

(27) 

c- A3 - 7x7 

IVaT-IKa; (28) 

De m§me que pour le cas de I'anche d mode unique, on va construire 
un modele num^rique dans lequel le deplacement de I'anche a Tinstant 
6chantillonne tn=:^est fonction de la pression acoustique, non pas ^ I'instant 

tn mais a I'instant tn.i= Cecl est possible du fait que la r6ponse 

impulsionnelle du modele est la somme de fonctions sinusoYdales amorties. 
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Pour remplir cette condition, on utilise une approximation de ioo sous la 
forme i© « f^Cz - 1) = f^{exp{\xn) - 1). 

De plus, pour ajouter au modele un parametre de contrdle 
suppl^mentaire, on remplace le coefficient B par un filtre tel que 

•3 » t_a exp(-iCT) ■ °" ajuster I'amortissement des harmoniques 

en fonction de I'amortissement de la fondamentale. Pour conserver la 
caracteristlque X(0)=1, les param^tres ba et ag sont lies par I'equation 

Le terme exp|^- i ^jest remplace par son equivalent ecliantilionn6 : 

z-°^=exp(-iOTDa) avec le retard Dg defini par : Da=E|^^j dans lequel E 
indique la partie entidre. 

1 

Si I'on note ^ P = ^ ^'^f ' ^"^^on transfert num§rlque s'6crit : 
X(z) _ C{zP^*' - a,z°" - b3z) 

Pe(z) (P + fe(z - ^)){z°^*' - a,z°^ + b^z) . ; 

^ C(z-^ -a^z-^ -b,z-°""') ' ■■ 



U - (fe(1 + - P)z-^ -ajp - fJz-2 + f.b^z-^^ - bjf, - p)z-°^-^ 
ce qui conduit a I'equation aux differences : 

X(n)=balPe(n-1)+ba2Pe(n-2)+baDlPe{n-Da-1) +aalX(n-1) +aa2X(n-2)+aaDX(n-Da)+aaDlX(n-Da-1) (29) 

dans laquelle les coefficients aai, aa2. aaoz, aaoi sont d6finls par : 

3 - fa(1 + a3)-p _ a,(p-fj b,(f, -p) 

aai , aa2 , ago—ba, aaDi=-^^-| — — 

et les coefficients bai, baz, baoi par : 
. C - Ca. . - Cb» 

bal= — , ba2= 2., baDI ^ 



fe fe " fe ■ 



On remarque que I'equation 29 ainsi §tablie permet de determiner le 
deplacement adimensionn6 x(n) de I'anche a I'instant echantillonne (n). a 
partir des Instants precedents. 

ii apparaTt done que le schema de calcul num6rlque du son dans le 
cas d'un modele d'anche multlmode peut dtre le mSme que celul de I'anche S 
mode unique en remplagant l'6quatlon 21 par I'dquatlon 29 qui est une autre 
transcription numerique du module (m) de la figure 4. 

A titre d'exemple, la figure 8, qui est analogue aux figures 1 et 6, 
donne, pour un module approche d'anche multimode calcute selon 
I'invention. deux diagrammes representant respectivement, en trait plein, d 
gauche la fonction de transfert et a droite ia reponse impulsionnelle avec une 
frequence d'echantillonnage fe=44100 Hertz, les parametres ayant des 
valeursfr=1 837.5 Hz. qr=0,2, aa=0. 

Sur les m§mes diagrammes, on a superpose, en traits Interrompus, la 
fonction de transfert et la reponse Impulsionnelle du modele de anche unique 
telles que representees sur la figure 6. 

Selon un autre developpement de I'invention, on peut aussi ameliorer 
le module de son de clarlnette de fagon a le rendre plus naturel en lui 
incorporant un certain bruit, le systdme §tant ainsi plus r§aliste. Etant donne 
que le bruit est cre6 par une turbulence au niveau de I'anche avant le d^but 
du tuyau, le bruit est ajoute d x(t). II apparaTt d'autre part que, en pratique, le 
niveau de bruit depend de la presslon du souffle alors que sa « couleur » 
depend de la pression des levres sur I'anche. En effet. d'un point de vue 
physique, plus I'anche est press^e. plus faible est I'ouverture entre I'anche et 
le tuyau et plus importante est la turbulence. On utilisera done un simple 
module de bruit dont le niveau est pllote par y et la brillance pilotee par ^. 

Pour la mise en oeuvre numerique d'un tel modele de bruit, on utilisera 
un fiitrage passe-bas d'un bruit blanc. La fonction de transfert de ce filtre est 
donn§e par T,(z) = MLz^^ |e coefficient bb 6tant pllote par y, et le 

coefficient ab pilote par ^. Les lois de variation de bb et aa peuvent etre 
detenninees de telle sorte que le son simule par le modele soit aussi realiste 
que possible. 

Les deux diagrammes de la figure 9 montrent a titre d'exemple, la 
variation du module du spectre de ia presslon acoustique externe 



correspondant au son prodult par !e module, respectivement sur le 
diagramme du haut pour une anche d mode unique et sur !e diagramme du 
bas pour une anche a mode multiple avec bruit additionnel, les param^tres 
de simulation 6tant les suivants : 

5 

fr=2205 Hz. qr=0.25, aa=0. y=0A4, ^=0.4. L=0.48 m, R=7.10-^m 

Comme on I'a Indique, le precede selon ['invention, tel qu'il vient d'etre 
decrit en detail, a pour objet la simulation de sons produits par un instrument 

10 de musique a anche et ^ resonateur cylindrique, du type ciarinette. Mais 
I'invention ne se limite pas a une telle application et peut, au contraire, faire 
I'objet de nombreux d6veloppements. 

En effet, a partir du modele physique non lineaire, faisant appel a des 
guides d'onde et sch6matis6 sur la figure 4, ainsi que de sa transcriptiqn 

IS numSrique selon la sequence des Equations 21 a 25, il va gtre possible,,-de 
simuler le fonctlonnement d'un resonateur ayant une geometric quelconqUe, 
en modifiant le modele d'imp6dance pour resonateur cylindrique et. les 
equations aux differences associ^es, pour les remplacer par des Imp^danoes 
et des equations aux differences plus complexes. ^ 

20 Tout en conservant les proprietes du modeie qui vient d'etre de<^t, 

d'autres modeles d'impedance plus complexes peuvent, en effet, etre b&as 
en combinant certains elements d'impedance en parallele ou en serie et en 
proposant des approximations numeriques permettant un usage explicite des 
variables physiques et un controie plus flexible de I'lnstrument numerique. . . 

25 A titre d'exemple, on va malntenant decrire certains developpements 

du modele de base pour resonateur cylindrique, en se referant aux figures 10 
a 13 qui representent des schemas de calcul equivalents faisant appel a des 
guides d'ondes et con-espondant ^ des resonateurs ayant des geometries 
diverses. D'une fagon generate, sur ces schemas, I'operateur C(©) 

30 represente I'lmpedance d'entree et C'\g>) I'admittance d'entree d'un 
resonateur cylindrique, le modeie numerique correspondant a C'\(o) etant 
obtenu en changeant simplement le signe du coefficient bo. 

Un premier perfectionnement du modele de base qui vient d'etre decrit 
en se referant aux figures 2 et 4, va permettre, en faisant appel, de fagon 

35 . analogue, a des guides d'onde, de realiser un modeie physique pour 
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rdsonateur cylindrique avec fmp§danoe termlnale. Un tel §l6ment permettra, 
par exemple, de relier entre elles des portions de r§sonateurs cylindrlques 
ayant des longueurs et des sections differentes, de fagon d simuler 
I'impedance d'entree d'un conduit de section variable, ou bien de prendre en 
compte I'impedance de rayonnement. 

Pour cela, on considere le formalisme de la ligne de transmission llant 
la pression et le d^bit acoustiques ^ I'entree du r^sonateur (Pe(<B), Ue(G))) et a 
son extremlte ouverte (PsCco), Us(q>)). 

DC 

Si I'on note = rimp6dance caracteristique on a : 

Pe(co) = cos(k(<3())L)Ps(oo) + iZcSin(k(G))L)Us((B) 
Ue((B) = ^ sin(k(ca)L)Ps({o) + cos(k(co)L)Us(co) 

En notant Z3 = I'imp6dance de sortie et R(a>) = ~ ^s(co) 

UsC©) + Z,((d) ' 

on peut ecrire rinnp6dance d'entree de deux fagons diff6rentes : 



Pe(03) Z, 



+itan(k((Q)L) 



C 



c e\ / 1 + 1 y ^ tan(k(o))L) 



(30) 



20 = 1 - R(a>) exp(-2ik(o))L) 

1 + R((D)exp(-2ik(a))L) 
L'equation 31 montre done que I'impedance d'un r6sonateur 
cylindrique avec impedance terminate peut §tre obtenue a partir de 
I'impedance d'un resonateur cylindrique sans impedance tenninale, en 
remplagant : 

25 exp(-2ik(cD)L) par R(co)exp(-2lk(co)L). 

La Figure 10 donne un schema equivalent de calcul faisant appel S 
des guides d'onde. pour la mise en ceuvre de l'equation 30, permettant de 
calculer I'impedance d'un resonateur cylindrique avec impedance terminate. 
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Un tel models permet de g6n6rer en cascade rimpedance d'entree 
d'un conduit ayant una geom6trie quelconque et peut m§me s'appliquer a la 
simulation du conduit vocal. 

Mais rinvention peut faire Tobjet de bien d'autres developpements et 
5 pennet en particulier, a partir du modele de base physique pour r6sonateur 
cylindrique schematise sur la figure 4, de construire des modeles specifiques 
pour la simulation de divers instruments de musique. 

C'est ainsi que dans un premier developpement du modele de base, 
on va construire un module pour resonateur conique, utilisable, par exemple, 
10 pour la simulation d'un saxophone. 

Si Ton appelle L la longueur du tuyau, R son rayon d'entree, 9 son 

ouverture, la distance x© entre le sommet du c6ne et Tentree est : 
R 
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Xe = 



(sin I] 



15 Ulmpedance d'entree, relative d rimp§dance caracteristique ...f 

= est alors dbnnee par Texpressidn : / 

TcR^ . 

- ' - (32) . 



j^x^ itan(k((a)L) 
c 

qui peut s'6crire sous ia forme : 



Pe(c«) 



1+ ^ 



itan(k(<»)L) 

La Figure 1 1 donne un schema equivalent de calcui faisant appel a 
des guides d'onde, dans lequel l'6l§ment not§ D est rop§rateur dlff6rentiel 
D=i<x) et I'element not§ C'\w) correspond au schema de la figure 2 en 
25 remplagant -exp (-2ik((o)L) par +exp (-2lk(oi))L). 

Las equations aux differences correspondantes peuvent etre Stablies 
en suivant un processus analogue a celui qui a §te decrit plus haut. 
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AinsI, en utilisant la transform§e billn^aire pour approximer io, la 
fonction de transfert numerique d'un resonateur conique est donn6e par : 

Pe(z) 1 



ZcUe(z) z + 1 1 - a,z~^ + bpZ 

2fe^(z-1)"'l-^i^"-ba2"" 
c 



1 1 
En notant : Gp = 1+ — et Gm = 1 . la fonction de transfert 

2f e ^ 2L ^ 

^ c * c 

se r6duit a : 

Pe(z) ^ 1 - (a^ + 1)z-^ + a,z-2 - boZ-2° + bpZ-^""^ 

ZcUe(z) Gp -(a,Gp + G,)z-^ +a,G^z-^ +boG„z-^ -boGpZ-^"-^ 

10 d'oD Ton tire I'equation aux differences : 

Pe(n) = bCoUe(n) +bciUe(n-1)+bc2Ue(n-2)+bcDUe(n-2D)+bCDiUe(n-2D-1) 
+ aciPe(n-1)+ac2Pe(n-2)+acDPe(n-2D)+acDiPe(n-2D-1) (33) 

15 dans laquelle ies coefficients bco, bci, bca, bcp et bcpi sont d^finis 

par: 

bco=J-. bc,= -^. bca=|l-, bco=-^. bcDi=^ 
°P Gp Gp 

20 et Ies coefficients aci. aca, aco et acoi sont d§flnis par : 

ac, ac =^ ac - ac -b 

Gp Gp 

De fegon analogue, I'lnvention peut s'appliquer au cas de r6sonateurs 
courts qui apparaissent. par exemple, dans I'embouchure d'un cuivre ou d'un 
25 trou de reglstre. 

Pour cela, on admet que le rayon du resonateur court est 
suffisamment grand pour conserver le module de perte utilise jusqu'a 
present. 



Cette approximation d'un resonateur court se fera en approximant 
['impedance Zi(cd)= i tan (k(co)l) pour des valeurs faibles de k(cD)l. 
On arrive aihsi S Texpression 

Z|((d)= i tan (k(o>)l) = G((d) + icDH(a)) (34) 

( 

. 1-exp 



-acj— 



dans laquelle G(cd) = — — . 

1 + exp 



V 



acJfl 



etH(Q))=-(l-G(a))), 



Dans un autre developpement du modele de base pour resonateur 
cylindrique, I'invention permet ^galement de simuler I'ensemble d'un . 
10 resonateur de type cuivre. 

Dans ce cas, on modelisera I'embouchure par un resonateur de 
Helmholtz comprenant une cavite hemispherique couplee avec un tuyau 
court cylindrique et un resonateur principal a tuyau conlque, 

Uimpedance d'entr§e de Tensemble du resonateur peut s'exprintpr 



15 par 



Ze(co) = 



Jl_ 



. V 1 

10) — - + 

pc^ . xe 



Zi tan(ki(G))Li) + Z 



2 



ICO 

c 



. xe 

ICO — 

1 + ^ 

itan(k2(<o)L2) 



4 ^ 

dans laquelle y- — nRl est le volume de la cavite hemispherique, Li 
o 

est la longueur du tuyau court, L2 est la longueur du tuyau conlque, Zi et Z2 
sont les impedances caracteristiques des deux tuyaux qui dependent de 
20 leurs rayons, ki(o)) et k2(o) tiennent compte des pertes et du rayon Ri et R2 
de chaque tuyau. 

Cette equation peut etre ecrite sous la forme : 
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Z„Z,((D) = 



iZ, tan(ki(©)Li) + Zj 



ICO 



ICO 



1 + 



itan(k2(co)L2) 



1 + ico 



V 



iZ^ tan(k^(o))Li) + Zj 



1(0 



100 



1 + 



itan(k2 (00)12) 



10 
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On peut ainsi etablir le schema equivalent de calcul represente sur ia 
figure 12, dans lequel Topirateur note Ci(cd) correspond au schema de la 
figure 4 et I'op^rateur note S(co) correspond d I'impedance d'entree du tuyau 
conique et au schema de la figure 11. 

Des m^thodes de resolution analogues a celles qui ont 6t§ decrites 
pr§cedemment permettent d'exprimer le module Equivalent num§rlque et les 
Equations aux differences conrespondantes. L'impedance du tuyau court peut 
etre modElisee en utilisant ['approximation exprimEe par I'Equatlon (34) 
precEdente. L'impedance du tuyau conique est representee par le modele 
correspondant ^ I'equation aux differences (33). L'admittance de la cavft6 

V 7-1 

ico — - est approximee par la transfonnation bilinEaire d- — oii d=2fe. 

pc^ z + 1 

En considErant I'association de la cavite hemispherique et du tuyau 
court comme un reso nateur de Helmholtz ayant pour frequence de 

resonance = c^^^ , on peut utiliser cette frequence coh pour approximer 

G(co) par G(o)h) et H(co) par H(ojh)- Les deux frequences utilisees pour le 

calcul des coefficients a^ and bo sont co, = c(127iL + 97c% +16L) 

' 4L(4L + 37cXe + 4x J ^ 
correspond au premier pic d'impEdance du tuyau conique, et coa = coh- De 

DC 

plus, on utilisera Z„ = = Z2 pour normaliser rimpedance d'entree. 

L'impedance num^rique d'un resonateur de type cuivre est ainsi 
donnee par I'expression : 
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s, 



n 



^« = Vd(z - 1) , 1 

^1 ^2 



qui peut etre simplifiee sous la forme : 

S^"^ he 7'^ + T^°^ hr z"^ 



Ze(z) = 



aCp — Zk«i 3C|^2 ^ — Skoo 3^Dk^ 



d'oD I'on tire Tequation aux differences : 

Pe(n)=5]bG,Ue(n - k) + 'fbCoKUeCn - k - 2D) (3f) 

k=0 k=*0 
k=4 k«3 

10 + EQ^KPe(n -k) + 2^CDkPe(n - k 2D) 

k«1 k=0 

?^ 

Les coefficients decoulent d'un calcul direct a partir de Tequation (35);r 

L'inventlon peut encore dtre appliqu§e ^ ia modelisatlon d'uh 

15 r6sonateur cylindrique avec trous de registre. On utilisera, a cet effet, d&k 

elements falsant appel a des guides d'ondes et correspondant 

respectivement a un meddle physique de tuyau cylindrique avec impedance 

terminate representant un tuyau de longueur L1 entre I'embouchure et le trou 

de registre, un modele de tuyau court qui represente le trou de registre de 

20 longueur hi et le rhodele de base pour tuyau cylindrique representant un 

tuyau de longueur L2 entre le trou de registre et Textremite ouverte. 

Uimpddance terminate de la premiere partie du tuyau peut s'6crire : 

1 



Zs(a))=- 



1 1 
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En ecrivant cette expression sous la forme : 
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1 , Za C,(a)) 
Z. C,((D) 

on peut etablir un schema general de calcul par guides d'onde 
represente sur la Figure 13. qui, dans le cas qui vient d'etre decrit, est utilis§ 
pour calculer I'lmpedance terminate par une combinaison en parallete des 
5 impedances ZaCzCco) et ZtCt(ffi). 

Dans le cas limite d'un trou de registre ferme, avec deux parties du 
tuyau ayant le meme rayon, on pose ZtCt(a))=oo et Z2=Zc 

d'oD il vient Zs(o>)=iZc tan(k(a))L2) et Zo(G))=iZc tan(k(co)(Li+L2)). 

10 L'impedance totale d'entree du tuyau peut alors s'exprimer par : 

Pe(«)) ^ (CaCcp) + C,(cD))ZtCt(Q) + Z,C,(a))Ca(a)) 
ZcUe(a)) (1 + C,(©)C2(co))Z.Ct((D) + Z^C^Co) ^ ^ 

En utilisant iffl=fe(1-Z"^) dans I'imp6dance du tuyau court, comme 
15 pr6c§demment, il en d^coule directement I'^quation aux differences. 

AinsI, le moddle de simulation du r^onateur cylindrique du type 
clarinette, obtenu par transposition directe des Equations slmpllfi^es du 
comportement physique de I'instrument, peut §tre adapts S la simulation 
d'instruments a resonateur non cylindrique, tel que le saxophone, la 
20 trompette ou autres Instruments a tuyau de cuivre. 

Mais I'invention ne se limite pas a un tel mode de realisation et aux 
adaptations qui viennent d'etre decrites car, sans s'ecarter du cadre de 
protection defini par les revendications, eile peut s'appliquer a la simulation 
d'autres types d'instruments, par exemple a corde frott§e comme le violon ou 
25 frapp^e comme le piano, les variables de la relation d'impedance etant alors 
la force exercee en un point sur la corde par I'archet ou le marteau et la 
Vitesse de la corde en ce point 

En effet, dans le cas d'une corde, la relation d'impedance mecanique 
entre la force exercee sur la corde en un point quelconque et la Vitesse en 
30 resultant, pemiet d'exprimer I'admittance mecanique de la corde en ce point 
sous forme d'une combinaison en paralldle des admittances de chaque 
portion de corde selon le schema de la figure 13, dans lequel les §l6ments 



notes Cz et Ct representent radmittance de chaque portion de corde et 
correspondent au schema de la figure 2, le filtre F(co) du modele de base 
pouvant etre exprime, soit a partir d'un modele de propagation d'ondes de 
flexion dans une corde munie de raideur, soit a partir de mesures 

5 experimentales. En effet, Tadmittance mecanique de chaque portion joue le 
role d'une impedance acoustique. 

En outre, la methode selon ['invention pourrait aussi s'appliquer a la 
simulation d'autres phenomenes oscillants, grace a une adaptation de 
certaines equations aux differences et un choix de coefficients tenant compte 

10 des caract§ristiques physiques des phenom6nes a simuler. 




REVENDICATIONS 

1. Proced6 de simulation num§rique d'une interaction non Iin6aire 
entre une source excltatrlce et une onde dans un resonateur, au moyen 
d'outils de calcul de signaux num^riques a partir d'equations dont la solution 

5 correspond a la manifestation physique d'un phenom^ne a slmuler qui se 
traduit, a chaque instant et en chaque point du resonateur, par une relation 
d'impedance entre deux variables representatives de I'effet et de la cause 
dudit plienomene a simuler. precede dans lequel on transcrit directement 
I'equation de I'imp6dance sous forme d'un modele faisant appel a des guides 
10 d'ondes pour realiser une interaction non lineaire entre les deux variables de 
la relation d'impedance. 

2. Proced6 selon la revendication 1, de simulation d'un phenomene 
oscillant, caract6ris§ par le fa'rt que le modele comporte, d'une part au moins 
une partie lineaire representant directement rimpSdance d'entr§e du 

15 resonateur et, d'autre part, une partie non lineaire mod§lisant le r6le de la 
source excitatrice du ph§nom§ne a simuler. 

3. Procede de simulation selon Tune des revendications 1 et 2, pour la 
synthase num6rlque, en temps r6el, d'un phenom6ne oscillant, caract6ris6 
par le fait que. d partir d'un systeme d'equations entre au moins deux 

20 variables representatives du comportement d'une onde dans le resonateur, 
on etabiit une expression de I'imp6dance d'entr6e du resonateur en tenne de 
guide d'onde lineaire, sans decomposition en onde aller-retour, de fagon a 
realiser au moins une partie lin§aire du modele. 

4. Procede selon la revendication 3, caracteris§ par le fait que la partie 
25 lineaire du modele est couplee a une partie non lineaire faisant intervenir 

revolution de la non-linearit§ telle qu'elle s'exprime entre les deux variables 
de la relation d'impedance d'entree du resonateur. 

5. Procede selon la revendication 4, caract6rise par le fait que la partie 
lineaire du modele de simulation numerique de I'equation de I'impedance fait 

30 appel S deux guides d'ondes 6l§mentalres r6alisant une fonction de transfert 
entre les deux variables de la relation d'impedance. 

6. Proc6d6 selon la revendication 5, caracterise par le fait que la partie 
lineaire a deux guides d'onde du module est coupl6e d une boucle reliant ia 
sortie a I'entree de ladite partie lin§aire et comprenant une fonction faisant 

35 intervenir la non-llnearite telle qu'elle s'exprime physiquement. 



7. Procede selon la revendication 6, caract§ris6 par le fait que le 
modele est pilote par au moins deux parametres representatifs de 
rinteraction physique non-lineaire entre la source et le resonateur, au moyen 
d'une boucle reliant la sortie a I'entree de la partie lineaire et comprenant une 

5 fonction non-lineaire jouant le role de source excitatrice pour le resonateur. 

8. Procede selon Tune des revendications precedentes, pour la 
synthese en temps reel de sons produits par un instrument a vent, 
caracteris§ par le fait que les deux variables de la relation d'impedance sent 
la pression et le debit acoustiques a Tentree du resonateur. 

10 9. Procede selon la revendication 8, caracterlse par le fait que. pour 

un resonateur cylindrique ayant une extr^mite ouverte, la partie linSaire du 
modele de transcription numerique de {'equation de i'imp6dance constitue la 
somme de deux guides d'ondes Slementaires ayant pour source excitatrice |e 
debit a Tentree du resonateur et realise la fonction de transfert : 
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Z (to) = = ^! exp(-2ik(co)L) 

^ Ue(co) 1 + exp(-2lk(o))L) . 1 + exp(~2ik(a))L) ^ 



dans laquelle : ^. 

- 00 est la pulsation de Tonde ^, 
20 - Ze{co) est Timpedance d'entree du resonateur, 

- Pe(cD) et Ue(ca) sont les transformees de Fourier des valeurs 
adimensionnees de la pression et du debit a Tentree du 
resonateur, 

- k(o>) est une fonction de la pulsation de Tonde qui depend du 
25 phenomdne a simuler, 

- L est la longueur du resonateur. 

10. Procede selon la revendication 9, caracterls6 par le fait que 
chacun des deux guides d'ondes fait intervenir un filtre ayant pour fonction 
30 de transfert : 

-.F(co)' = -exp(-2ik(a>)L) 



et repr§sentant un trajet aller-retour d'une onde avec changement de 
signe d I'extr^mite ouverte du r6sonateur, chaque guide d'onde 
correspondant S un terme de ('Equation de rimp^dance. 

11. Precede selon la revendication 10, caracterise par !e fait que le 
module est pilote par la longueur du resonateur et au moins deux paramdtres 
representatlfs de I'lnteraction physique non lineaire entre la pression et le 
debit ^ I'entree du resonateur, au moyen d'une boucle reliant la sortie ^ 
I'entree de la partie lineaire et comprenant une fonction non lineaire jouant le 
role de source excitatrice pour le resonateur. 

12. Precede selon la revendication 11, caracterise par le fait que la 
fonction non lineaire a comme parametres d'entree la pression et le 
deplacement de I'organe de formation des vibrations et est contrdlee par au 
moins deux parametres de simulation du jeu d'un instrumentiste. 

13. Precede selon la revendication 12, caracterise par le fait que les 
parametres de jeu pour le contr6le de la fonction non lineaire sont : 

- un paramdtre <; caracteristique de {'embouchure et de Taction de 
rinstrumentiste sur I'organe de formation des vibrations, 

- un paramdtre y repr^sentatif de la pression appliqu6e sur I'organe 
de fonnatlon des vibrations. 

14. Proc6d6 selon I'une des revendications 10 a 13, caracterise par le 
fait que, pour la synthase en temps reel des sons a simuler, on realise une 
formulation dans le domaine temporel de la rSponse Impulslonnelle de 
I'impedance du resonateur, en approximant les pertes representees par le 
filtre au moyen d'un filtre numerique approche. 

15. Precede selon la revendication 14, caracterise par le fait que I'on 
utilise un filtre numerique a un pole en exprimant la fonction de transfert 
approchee sous la forme : 



1 - a, exp(-iiir) 



(13) 



dans laquelle : 

■ ^ = 7^ . /e etant la frequence d'echantillonnage, 



e 
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• D - X — est le retard pur correspondant ^ un trajet alter ou retour 

c 

de Tonde dans le resonateur, 

les coefficients bo et ai sont exprimes en fonction des 

F(<jo)^ = W{js)\ pour une vateur 

5 de la pulsation correspondant a la frequence fondamentale de jeu et une 
autre valeur ©a correspondant a un harmonique. 

16. Precede selon la revendication 15, caract6rise par le fait que les 
coefficients bo et ai sont obtenus par resolution du systeme d'equation 



10 I (1 + ^1 - 2a^ cos(tn^)) = b; 

|F(©,)2|'(1 + a? - 2a, cosCtrr^)) = b^ . T 

avec |F((jo)^j^ = exp(-2aC-^L) , lesdits coefficients etant donnes p^r 



15 lesformules 



A1 - A2 - V(A, - A,)' - (F, 



bo = 



^2F,F,(c, - c,)(A, - A, - ^A, - A,f - (F, - F,)^) 

F, -F2 

dans lesqueiles 

C, = COS(tEfi),C2 = C0S{OT2).Fi = |F((B1)^|^F2 = |F((D^)^|^ 

A, = FjC,, = f"2^2 

17. Precede selon la revendication 16, de simulation d'un instrument a 
25 anche a resonateur cylindrique, d partir d'une modelisation physique regie 
par le systeme d'equatlon : 



1 



x(t) q, dx(t) 
It^ dt 



+ x(t) = p,(t) 



Pe(a>) 



i tan 




5 



Ue(t) 



1(1 - sign(Y - x(t) - 1))sign(Y - Pe(t)K(1 - Y + x(t))^|y - pjt)| 



dans lequel ©r est la frequence de resonance et qr est le facteur de qualite de 
I'anche, caracterise par le fait que ledit systdme d'equations est r§solu dans 
le domaine temporel ^ partir d'une formulation echantillonnee equivalente de 
la r^ponse impulslonnelle du deplacement de I'anche et de la relation 
10 d'imp6dance qui se tradult par le systeme d'equations : 

x(n)=biaPe(n-1)+aiaX(n-1)+a2aX(n-2) (1 8) 

Pe(n)=Ue(n)-ai uq (n-1 )-boUe(n-2D)+aiPe(n-1 )boPe(n-2D) (1 9) 



lesdltes equations etant utilisees de maniere s6quentielle en 
regroupant les termes ne dependant pas de I'echantillon temporel n. de 
fapon a calculer successivement : 



15 



"e(n) = s'g"(Y - x(n) - 1))sign(Y - Pe(n)Kn)- y + x(n))^|Y - p, (n)| (20) 



x(n)=biape(n-1 )+aiaX(n-1 )+a2aX(n-2) 



(21) 



V=-aiUe(n-1)-boUe(n-2D)+aiPe(n-1)-boPe(n-2D) 



(22) 




(23) 



Ue(n) = 2 signCY - V)(-bCoW2 + \N^{bc^\Nr + 4|y - V| 



(24) 



Pe(n)=boCoUe(n)+V 



(25) 



18. Proc6d§ selon la revendlcation 17, pour une simulation plus 
30 r^ailste du son produit, caracterise par le fait que, en negligeant le 



rayonnement, la pressfon externe s*exprime comme la d6rivee en fonction du 
temps, du debit sortant, sous la forme 

Pext(t)=-|(Pe(t) + Ue(t)) (26) 

5 et est calculee, ^ chaque instant echantilionne (n), par difference entre 

les sommes de la pression interne Pe et du debit Ue, respectivement ,d 

Tinstarit (n) et a Tinstant (n-1). 

19. Procede selon la revendication 17 de simulation d'un instrument a 

anche multimode, caracterlse par le fait que le calcul de la pression 
10 acoustique et du debit au niveau de i'embouchure s'effectue par resolution 

sequentielie d'un systems d'equations dans lequel le deplacement de Tanche 

a chaque instant (n) est de la forme : 



15 



20 



x(n)=bai Pe(n-1 )+ba2Pe(n-2)+baoi Pe(n-Da-1 ) +aaix{n-1 ) +aa2X(n-2)+aaDX(n-Da)+aaDiX(n-Da-1 ) . ' 
les coefficients aai, aa2, eaD2j 3aD1 etant definis par : 



«a1 Z 1 aa2 ~ % aaD — Oa, oaDI 7 

et les coefficients bai , ba2, baoi par : 



bal=: ^ , ba2 =: — 1 baDI 



a 



fe' fe fe 

2 1 

en posant : = . . A2=— co^q^ et A3=A2Aiqr. 

O^rV^r +4 2 

25 les equations suivantes etant les memes que pour une anche d mode 

unique. 

20. Procede selon Tune des revendications 17 a 19, caracterise par le 
fait que Ton realise un models pour resonateur cylindrique avec impedance 
terminate a partir du modele de base correspondant a un resonateur 
30 cylindrique et constituant la somme de deux guides d'onde faisant intervenir 



chacun un filtre ayant pour fonction de transfert -F(o)^=-exp(-2ik(o))L), en 
rempla9ant ('expression exp(-2ik(o>)L) par I'expression R(cD)exp(-2ik(a))L), 

dans laquelle R(co) = ~ ^"^^^ 

Zc +23(0) 

Zc 6tant rimp§dance caract6ristique etZg rimp6dance de sortie .^1^ 

5 21. Precede seion Tune des revendications 17 a 20, caracterise par le 

fait que, a partir du modele d'impedance pour un instrument a resonateur 
cylindrique et des equations aux differences associees, on construit d'autres 
modeles d'impedance plus complexes pour la simulation de phenom^nes 
osdllants produits par un r6sonateur de fonne quelconque en combinant des 

10 Elements d'impedance en parallele ou en serie et en utilisant des 
approximations numeriques pour un usage explicite des variables physiques 
intervenant dans la production desdits ph6nom§nes oscillants et un contrSIe 
plus flexible du resultat de la simulation. 

22, Precede selon Tune des revendications 17 a 21, caracterise par le 

15 fait que, d partir du modele de base pour resonateur cylindrique dans lequel 
la reponse impulsionnelle du deplacement de Tanche se traduit par un 
systeme d'equations aux differences donnant, a chaque instant (n), le 
deplacement x(n) la pression pe(n) et le debit Ue(n) a Tentree du resonateur, 
on construit un modele pour resonateur conique dans lequel r§quation de la 

20 pression est de la forme : 

Pe(n) = bCoUe(n) +bciUe{n-1)+bC2Ue(n-2)+bcDUe(n-2D)+bCDiUe(n-2D-1) 

+ aciPe(n-1)+ac2Pe(n-2)+acDPe(n-2D)+acDiPe{n-2D-1) (33) 



dans laquelle les coefficients bco, bci, bc2, bco et bcpi sont definis 



25 par : 



bCo=-^. bCi=-ii±l. bC2=-^. bCD=-^. bCD1=-^ 

C5p Gp Gp Gp Gp 

et les coefficients aci, ac2, aco et acoi sont definis par : 

30 ac _ QjQp Qm _flGm _ ^O^^m _u 

ju - — , ac2 - -— , acp - — — , ac^, - 

Gp Gp Gp 
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1 1 

en notant : Gp = 1+ et Gm = 1- 



2f , ^ 2fe ^ 

23. Proc6de selon Tune des revendications 17 a 21, caracterise par le 
fait que. a partir du modele de base pour resonateur cylindrique, on construit 
un module pour un resonateur court ayant une longueur I, en approximant 
5 rimpedance selon Texpression : 

Zi(g>)= i tan (k(a))l) s G(©) + mH(<o) (34) 



1 - exp 

dans laquelle G(cd) = . 

1 + exp 



- ccc J| I 



etH(o))=-(l-G(G>)). 
c 



10 24. Proc6d6 selon rensemble. des revendications 21, 22, 23, pourila 

simulation d'un instrument du type culvre, caracterise par le fait que/pn 
model Ise {'embouchure par un resonateur de Helmholtz comprenant iipe 
cavit6 h6mispherlque couplee avec un tuyau court cylindrique et^^'un 
resonateur principal a tuyau conique, I'impedance d'entree de I'ensemble^Ju 

15 resonateur pouvant s'exprimer par I'expression : ? 

_1_ 

Z„ 



Ze(ca) 



. V 1 

10) - + 



pc^ . X. 

^ !<0 — 



tan(ki(<D)L,) + Zj 



Xei 

ICO — 



1 + 



itan(l<2(o>)L2) 



dans laquelle V=-^7tRb est le volume de la cavite h6mispherique, Li 
6 

est la longueur du tuyau court, L2 est la longueur du tuyau conique, Zi et Z2 
20 sont les impedances caract6ristiques des deux tuyaux qui dependent de 
leurs rayons, ki(o)) et kaCoo) tiennent compte des pertes et du rayon Ri et R2 
de ciiaque tuyau, et que. d partir du modele de base pour resonateur 
cylindrique, de ses extensions au tuyau conique et au tuyau court, on 
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construit un moddle pour resonateur de type cuivre en exprimant la pression 
a I'embouchure par I'equation aux differences : 



25. Proc6d§ selon la revendlcation 21, pour la simulation d'un 
ph6nom§ne oscillant dans lequel ies deux variables physiques de la relation 
d'lmp^dance sont la force exerc6e en un point d'un systeme mecanique tel 
qu'une corde engendrant des vibrations et la vitesse en ce point, caracteris6 
par le fait que Ton exprime I'admittance en ce point sous la fonne d'une 
combinalson en parall6le des admittances de chaque portion de corde, de 
part et d'autre dudit point, chaque admittance mecanique 6tant obtenue d 
partir du modele de base d§crivant rimp§dance acoustique d'un resonateur a 
tuyau cylindrlque, en exprimant la vitesse au point considere de la corde en 
fonctlon de la force exercee en ce point, le filtre F(q)) du module de base 
pouvant etre exprim§, soit a partir d'un module de propagation d'ondes de. 
flexion dans une corde munie de raideur, soit a partir de mesures 
experimentales. 



I<=4 k-3 



Pe(n)=2]bc,u^(n - k) + ^bCokUeCn - k - 2D) 



(36) 




+ EaCkPe(n - k) + XacDKPeCn - k - 2D) 
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